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Abstract 

Spectral studies of a luminous X-ray source, NGC 253 X21, are presented based on archival Chandra and 
XMM-Newton data. The Chandra observation on 1999 December 16 detected the source at a bolometric 
luminosity of 0.3 x 10 39 erg s _1 (assuming isotropic emission), while an XMM-Newton observation on 2000 
June 3 revealed a short-term source variation in the range of (0.6-1.3) x 10 39 erg s _1 . All spectra from these 
observations were successfully modeled by emission from an optically thick accretion disk. The average 
inner disk radius was kept constant at 63- (cos 60°/ cos i) 1 / 2 km, where i is the disk inclination, and did 
not vary significantly, while the disk inner temperature changed in the range of 0.9-1.4 keV. Assuming 
that this object is an accreting Schwarzschild black hole, and that the disk inner radius coincides with its 
last stable orbit, the mass of the black hole is estimated to be ~ 7 Mq. The disk luminosity corresponds 
to (30-120)-(cos 60°/cos i) % of the Eddington limit of this black hole. Therefore, this luminous X-ray 
source, NGC 253 X21, is understood consistently to be an accreting stellar mass black hole in the standard 
disk state. 
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1. Introduction 

Since the era of the Einstein Observatory, many com- 
pact X-ray sources have been detected in nearby galax- 
ies (e.g., Fabbiano 1989). Sources with X-ray luminosi- 
ties below ~ 2 x 10 38 erg s^ 1 , which corresponds to the 
Eddington limit for a neutron star of 1.4 Mq, can be ac- 
creting neutron stars, while those with higher luminosities 
can be considered as black hole binaries (hereafter, BHB). 
The most luminous class of compact objects, known as 
ultra-luminous compact X-ray sources (ULXs; Makishima 
et al. 2000), show X-ray luminosity of ~ 10 39 -10 40 erg s _1 . 
Supposing that ULXs obey the Eddington limit, their high 
luminosities imply that they are accreting black holes as 
massive as 10— 100 Mq or more. Such luminous X-ray 
objects, which are rather rare in our Galaxy and LMC, 
are suggested to reside preferentially in galaxies with high 
star-formation activities (Gilfanov et al. 2004). Therefore, 
starburst galaxies provide a good opportunity to study 
stellar black holes of various masses. 

In studying Galactic/Magellanic BHBs, intensity- 
correlated spectral changes have been used as a key di- 
agnostic tool. When the luminosity is above ~ 3% of the 
Eddington limit, a BHB is often in the standard-disk (al- 
ternatively called high or soft) state, wherein the X-ray 
spectrum is successfully described by a so-called multi- 

* Partly based on the observations obtained with XMM-Newton, 
an ESA science mission with instruments and contributions di- 
rectly funded by ESA Member States and NASA, USA. 



color disk model (MCD model; Mitsuda et al. 1984), pre- 
dicted by the theory of optically thick standard accretion 
disks (Shakura, Sunyaev 1973). The bolometric disk lu- 
minosity, Ldisk, varies as oc Tj n 4 , where T- m is the highest 
color temperature of the disk. Equivalently, the calcu- 
lated innermost disk radius, i?i n , stays constant, at a value 
consistent with the radius of the last stable orbit around 
the black hole. As the source becomes more luminous 
(close to the Eddington limit), R ln apparently becomes 
variable, as R\ n oc Tj n _1 , probably because the disk makes 
a transition into a slim disk (Kubota et al. 2001; Kubota, 
Makishima 2004). A very similar behavior has been ob- 
served for several of the most luminous ULXs, which 
show MCD-like spectra (Mizuno et al. 2001; Sugiho 2003; 
Kubota et al. 2002). It is therefore of keen interest to look 
for such intensity-correlated spectral changes among ex- 
tragalactic BHB candidates, in particular including those 
with lower luminosities, and to examine whether the be- 
havior indicates standard disks or slim-disks. 

With this perspective, we analyzed the XMM-Newton 
and Chandra data of NGC 253 X21 (Vogler, Pietsch 1999; 
also named Source 5 in Fabbiano, Trinchieri 1984). The 
nearby Sc-type edge-on galaxy NGC 253 (at a distance 
of 2.58 Mpc; Puche, Carignan 1988) is one of the well- 
known starburst galaxies with prominent diffuse X-ray 
emission (e.g., Strickland et al. 2002). Many point-like 
X-ray sources have also been detected since the era of 
the Einstein Observatory. Among them, X21 is one of 
the brightest sources, located approximately 4' southwest 
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Fig. 1. The 0.5-10 keV light curve of X21 (in counts s" 1 ), 
averaged over the two MOS cameras of XMM-Newton. The 
background is not subtracted. 

from the center of this galaxy. Voglcr, Pietsch (1999) 
showed that this source is time-variable by analyzing 
ROSAT data, and Pietsch et al. (2001) analyzed the 
XMM-Newton data of X21 and reported that its lumi- 
nosity varied by a factor of ~ 2 during a 40 ks observa- 
tion. In a large sample of ULXs observed with Chandra 
and XMM-Newton (Sugiho 2003), this source showed a 
large luminosity variation in the range of ~ 10 38 erg s _1 - 
10 39 erg s _1 . Therefore, this object deserves a detailed 
study in comparison with Galactic/Magellanic BHBs. 

2. Observations and Data Reduction 



The XMM-Newton (Jansen et al. 2001) observation of 
NGC 253, to be utilized in the present paper, was carried 
out on 2000 June 3. The data, retrieved from the XMM- 
Newton archive, were processed using the Science Analysis 
System (SAS), version 5.4.1. To create response matrices, 
the software package rmfgen was used. Since X21 was 
located on the CCD gap of the PN detector (Striider et 
al. 2001) during the observation, the present paper utilizes 
only the MOS (MOS 1 and MOS 2; Turner et al. 2001) 
data. For each EPIC camera, events with pattern to 12 
were accumulated over a region of 15", centered on the 
X21 image. A total of 35 ks of good data were obtained 
for the two MOS cameras by discarding the time region 
of a short background time flare. Figure 1 shows a 0.5- 
10 keV light curve with an average count rate of each MOS 
camera of 0.06 cts s" 1 , which roughly gives a 0.5-10 keV 
flux of 1 x 10~ 12 erg s _1 cm~ 2 . As shown by the dashed 
line at the bottom, the background counts were negligible 
(~ 1% of the source counts in this range). The signal 
count rate varied by a factor of two, as already reported 
by Pietsch et al. (2001). 

The present work also utilized a Chandra (Wcisskopf 
et al. 2000) observation of NGC 253, performed on 1999 
December 16. The Chandra data consisted of a net ex- 



> 

M 10 -2 



CO 
CO 

o 
o 



10 



10 



-3 _ 



-4 



-2 



: r 

: -MOS 1 
- "MOS 2 


1 1 1 1 1 1 1 II 1 1 • • • ■ 1 • • ■ • 1 • ■ ■ ■ 1 • • l± 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 lllll 


mi nil nil nil nil — i — 


1 1 1 1 II 1 1 1 1 1 1 1 1 1 1 

1 


1 1 1 1 1 1 II Mil 

1 



0.5 0.7 1.0 2.0 3.0 5.0 7.0 
Energy [keV] 

Fig. 2. Time-averaged XMM-Newton spectra of X21, to- 
gether with the data residuals from the best-fit MCD model. 

posure of 14 ks, when the target source was on the S3 
chip. The data reduction was carried out using CIAO 
software version 2.3, and the response matrix was made 
with mkrmf. The source events were extracted from 
a circular region of 3". 7 radius around the source cen- 
troid. The average 0.5-10 keV count rate was 0.06 cts s _1 , 
which corresponded to a 0.5-10 keV energy flux of 3 x 
10 -13 erg s _1 cm~ 2 , a factor of three lower than the av- 
erage flux observed with XMM-Newton. There was no 
significant time variation beyond the Poisson noise. 

3. Analyses and Results 

3.1. Averaged Spectra 

Figure 2 shows the time-averaged spectra of X21 ob- 
tained with the XMM-Newton MOS cameras. The back- 
ground events were extracted from a source-free circu- 
lar region of 2' radius on the same chip as the source. 
Following the previous studies on ULXs (e.g. Makishima 
et al. 2000), the MCD model and a single power law (PL) 
model were fit to the spectra. As summarized in table 1, 
the MCD model successfully describes the observed spec- 
tra with jv = 149.9/152, while the PL model is less 
successful with \ 2 l v = 175.9/152. The value of the hy- 
drogen column density, Ah, obtained with the MCD fit, 
is roughly consistent with the line-of-sight Galactic value 
of A H = 1.4 x 10 20 cm" 2 (Dickey, Lockman 1990), while 
that with the PL model is by a factor of 20 higher than 
the Galactic value. Furthermore, the value of Ah esti- 
mated from the MCD fit is consistent with the value from 
a spectral analysis of diffuse emission in this region of 
NGC 253 reported by Strickland et al. (2000). Thus, the 
MCD model explains the data more appropriately than 
the PL model. Employing a correction factor for the inner- 
boundary condition of £ = 0.41 (Kubota et al. 1998) and a 
color-to-effective temperature correction factor of k ~ 1.7 
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Table 1. The best-fit parameters of the energy spectra of NGC 253 X21* 
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Fig. 3. Chandra ACIS spectrum of X21, together with the 
best-fit MCD model. The bottom panel shows the fit residu- 
als. 



(Shimura, Takahara 1995), the best-fit MCD model yields 
the innermost disk radius, i?i n = (59 ± 6) ■ Q 1 ^ 2 km with 
£ = cos 60°/ cos i, where i is the inclination of the system. 
By using a formula for the disk bolometric luminosity, 
Tdisk, by Makishima et al. (2000), 



Jdisk 



: 47T 



T 

-L 1 I 
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where a is the Stefan-Boltzmann constant, Tdisk is then 
calculated as 8.6 x 10 38 • C erg s _1 , which corresponds to 
the Eddington limit of a 5.8 M object. An MCD plus 
PL model with a fixed photon index of T = 2.0 gives essen- 
tially the same result within 90% errors, because the PL 
component in this two-component fit is not significant. 

The Chandra ACIS spectrum of X21 was analyzed in 
the same manner. Figure 3 shows the obtained source 
spectrum, after subtracting the background spectrum, 
which was extracted from a source-free circular region of 
30" radius on the same chip as the source. The spec- 



Fig. 4. Relation between the highest color temperature of 
the disk (T; n ) and the calculated innermost disk radius (ilj n ) 
of NGC 253 X21. The dashed line represents R in = 63 km. A 
disk inclination of i = 60° is assumed. 

trum was fit with the MCD model and the PL model, 
as was done for the XMM-Newton data. As shown in 
table 1, the spectrum can be described by either model 
because of the poor statistics. However, the MCD fit is 
again more reasonable from the viewpoint of Nu- The 
Chandra data thus yield i?i n = (69^^) ■ C 1 ^ 2 km and 
Tdisk = 2.9 x 10 38 ■ C erg s" 1 . Between the XMM-Newton 
and the Chandra observations, the value of i?i n was kept 
almost constant at ~ 60 • C 1 / 2 km, while Ldisk changed by 
a factor of three. 

3.2. Intensity- Correlated Change of the Spectrum 

As shown in figure 1 and reported by Pietsch et 
al. (2001), X21 exhibited a significant intensity variation 
by a factor of two during the XMM-Newton observation. 
Accordingly, the entire observational span was split into 
three periods; covering low-flux, middle-flux, and high- 
flux periods, as denoted by the three arrows in the figure. 
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Fig. 5. Scatter plot between the disk bolometric luminosity 
(Z/disk) and the highest color temperature (Ti n ) obtained from 
the MCD-modcl fit. All data points of NGC 253 X21 refer 
to the present work. Those of other ULXs, IC 342 Source 1, 
IC 342 Source 2, M 81 X-6 and NGC 1313 Source B, were ob- 
tained with ASCA (Mizuno et al. 2001). The data points of 
the Galactic BHB, XTE J 1550 - 564, are taken from Kubota, 
Makishima (2004), and are shifted downward for clarity by 
a factor of two. The inclination is assumed to be 60° for 
all data. The dashed lines represent Ldisk °c Ti n 4 , hence 
mean constant value of R ln , and the dot-dashed lines rep- 
resent I/disk 0(1 Tin 2 - The parameter rj specifies the disk bolo- 
metric luminosity normalized to the Eddington luminosity 
(rj = ^disk/TEdd)- F° r NGC 253 X21, the same symbols with 
figure 4 are used. 

Spectra constructed from each period were investigated 
with the absorbed MCD model to clarify the spectral evo- 
lution in response to the intensity change. As shown in 
table 1, this model reproduced each spectrum very well. 

Figure 4 shows the values of i?j n against T[ n , based on 
the intensity-sorted XMM-Newton spectra (filled inverted 
triangles) and the Chandra spectrum (filled square). It 
clearly shows that the value of i?j n remained constant dur- 
ing the significant variation of T- m . The data points are 
consistent with the behavior of BHBs in the standard-disk 
state, where R m remained constant. The average value of 
i?i„ is determined as (63+ 9 ) • C, 1 / 2 km. 

In figure 5, the calculated values of Ldisk are plotted 
against T- m for each period, in order to examine the spec- 
tral behavior against the luminosity. This diagram shows 
that the data points satisfy the relation Ldisk ^ ^fa 4 a s 
Ldisk changes over (0.3-1.3) xIO 39 ■£ erg s _1 . In fact, a 
logarithmic slope of 3.6±0.6 is obtained when these data 
points are fitted with a power law. 

4. Discussion 

Through the combined use of archival Chandra and 
XMM-Newton data, an overall luminosity change by a 
factor of three was observed from NGC 253 X21. The 
X-ray spectra at all intensity levels were successfully de- 
scribed by the MCD model, and the highest disk luminos- 
ity reached 1.3 x 10 39 • C erg s" 1 . This luminosity, on one 



hand, is the lowest among those of ULXs. On the other 
hand, this source is shinning at a higher luminosity than 
most Galactic/Magellanic sources. Except for a few exam- 
ples, such as GRS 1915 + 105 and XTE J1550 - 564 in its 
outburst peak, no compact X-ray source in our Galaxy has 
so far exhibited a luminosity higher than 1 x 10 39 erg s _1 
(e.g., Gierlihski, Done 2004; Zycki et al. 1999). 

In order to study the spectral behavior of X21 compared 
with luminous ULXs and XTE J1550 — 564, which is a 
representative Galactic BHB, their Ti n -Ldisk relations are 
plotted together in figure 5. While the luminosities of 
luminous ULXs are observed to vary as oc Tin , that of 
NGC 253 X21 varied as Ldisk oc T- m A . The latter relation 
is considered to be a signature of standard accretion disks 
(Shakura, Sunyaev 1973), as has been observed from a fair 
number of BHBs, including XTE J1550- 564, LMC X-3 
(Ebisawa et al. 1993), GS 2000 + 25 (Takizawa 1991), and 
GS 1124- 684 (Ogawa 1992). The present results thus 
provide one of the first confirmations of the standard-disk 
property in BHB candidates outside the local group of 
galaxies. 

As shown in Section 3, the value of R m of X21 was 
estimated to be 63 • C, 1 !" 1 km. The innermost stable or- 
bit of a non- rotating BH is given as 3R S , where R s = 
2GM/c 2 is the Schwarzschild radius. Thus, the mass 
of the central black hole of this source is estimated to 
be 7 • C 1 ^ 2 Mq. This is a very common value found 
among BHBs in our Galaxy and LMC (e.g., McClintock, 
Remillard 2005). Thus, the observed disk luminosity of 
(0.3-1.3)xl0 39 • C erg s" 1 corresonds to (30-120) • C % of 
the Eddington limit for the inferred BH mass. Therefore, 
NGC 253 X21 can be understood consistently as an ac- 
creting BH with an ordinary stellar mass, in which a stan- 
dard accretion disk is radiating at a fair fraction of the 
Eddington limit. In particular, the observed values of 
Ti n = 0.9-1.4 keV are reasonable for the inferred BH mass 
and luminosity, thus making X21 free from the problem 
of "too high a disk temperature" observed from the most 
luminous ULXs (Makishima ct al. 2000). 

According to the present results, the least-luminous 
class of the entire ULXs population may contain ordinary 
BHBs. Similar examples include M 33 X-8 (Takano et 
al. 1994; Makishima et al. 2000), NGC 2403 X-3 (Kotoku 
et al. 2000), and NGC 4449 Source 2 (Miyawaki et 
al. 2004). Future studies of luminosity-correlated spectral 
evolution with a larger sample of extragalactic X-ray 
point sources will elucidate the detailed composition of 
ULXs. 

The authors would like to thank Chris Done for her 
valuable comments. T. Tanaka is supported by research 
fellowships of the Japan Society for the Promotion of 
Science for Young Scientists. 
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